2023 Session B 


B1. Consider an m-by-n grid of unit squares, indexed by (i,j) with 1 <i < mand1 <j <n. 
There are (m—1)(n—1) coins, which are initially placed in the squares (i, j) with 1 < i < m—1 
and 1 < j <n-—1. Ifa coin occupies the square (i, j) with i < m — 1 and j < n — 1 and the 
squares (i+1,7), (i,j +1), and (¢+1,7+1) are unoccupied, then a legal move is to slide the 
coin from (i,j) to (¢+1,j +1). How many distinct configurations of coins can be reached 
starting from the initial configuration by a (possibly empty) sequence of legal moves? 


—2 
Answer: & ila ) 
m— 1 


Solution 1: Think of (1,n) as the northwest corner of the grid and (m,1) as the southeast 

corner. Consider the unoccupied squares. Initially, they consist of all the squares in the 

north row (j = n) and/or the east column (i = m). We think of these squares as forming 

a lattice path from (1,n) to (m,1), and we denote this path by E... ES... S, representing 
m—l1 n—l1 

the sequence of eastward and southward steps that traverse the path from northwest corner 

to southeast corner. 


Claim. The unoccupied squares always form a lattice path from (1,n) to (m,1), with a total 
of m—1 eastward steps and n — 1 southward steps, and all such lattice paths can be achieved 
by a sequence of legal moves. 


Proof. The squares (1,n) and (m, 1) always remain unoccupied, because no legal move can 
slide a coin to either square. A legal move changes an ES portion (consisting of three 
unoccupied squares that make the move legal) of an unoccupied lattice path to SE, so by 
induction, after every move the unoccupied squares continue to form a lattice path consisting 
of the same number of E and S steps. 

On the other hand, given a lattice path of E and S' steps from (1,7) to (m, 1), which must 
consist of m + n — 1 squares, consider the configuration with coins on every square not on 
the path. If the sequence of steps does not contain an SE portion, then we are in the initial 
configuration. Otherwise, choose any SE portion; the three squares connected by the SE 
steps must have coordinates (i, j +1), (i,j), and (i+ 1, j), and there must be a coin in square 
(i+1,j +1). Sliding that coin to square (i, j) is the reverse of a legal move, and changes SE 
to ES in the path. Continue to change instances of SE to ES in this way until the initial 
configuration is reached. Reversing the sequence of coin slides made yields a sequence of legal 
moves from the initial configuration to the given configuration. 


Thus, the total number of configurations is the number of possible lattice paths, which is 


the number of different sequences of m — 1 E’s and n — 1 S’s; this number is Guana 
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B2. For each positive integer n, let k(n) be the number of ones in the binary representation 
of 2023 - n. What is the minimum value of k(n)? 


Answer: 3 


Solution 1: Clearly k(n) must be more than one, since no power of 2 is a multiple of 
2023 = 7-177. Furthermore, powers of 2 are all congruent to 1,2, or 4 (mod 7), and no sum 
of two of these residues can be zero. So k(n) must be at least three. We now show that this 
can be achieved. 

Note that the binary representations of both prime power divisors of 2023 have exactly 
three ones, as 7 = 2?+2!+2° and 177 = 289 = 28+2°+2°. It is now possible to piece together 
three powers of 2 whose sum is simultaneously a multiple of 7 and 289 using the Chinese 
Remainder Theorem, since the multiplicative order of 2 modulo 7 is 3, which is coprime to 
16-17 = 272, the order of the group of multiplication modulo 289 (from which it follows that 
2°72 = 1 (mod 289), by Euler/Lagrange’s Theorem). 

In particular, the triple (0, 277, 8) is congruent to (0, 1,2) (mod 3) and (0,5,8) (mod 272), 
so 2277 4 28 + 2° is a multiple of 2023. (Another natural solution here is 278° + 2° + 2°.) 


Solution 2: As in Solution 1, k(n) > 3 for all n. This solution gives an alternative simple 
approach for constructing n that achieve k(n) = 3. 

Since 24 = 17 — 1, raising both sides to the 17th power and using the binomial theorem 
yields 268 = —1 (mod 172). Also, squaring both sides yields 2136 = 1 (mod 177). Since 
268 — 967 1 967 we therefore have 2°” 42°" +2° = 0 (mod 177). As in Solution 1, 2'+2/+2* = 0 
(mod 7) if and only if {i,j,k} = {0,1,2} (mod 3). Since 136 = 1 (mod 3), this can be 
achieved by 2203 + 267 + 20, which is a multiple of 2023. 

As a minor variant, we similarly have 26° + 2° + 29 = —2 + 1 + 1 = 0 (mod 172), which 
naturally leads to the following solutions: 29 +2196 +2272, 2205 4.9272.4.99 and 2341 4.9136 4.90. 


Solution 3: This solution provides an exhaustive description of all possible triples of non- 
negative integers (i, j,k) such that 2° +21 +2% = 0 (mod 2023). All solutions can be reduced 
to a canonical form by factoring out (and removing) common powers of 2, and without loss of 
generality the terms may be ordered as 2'+2/+2°, where i > j. The Chinese Remainder The- 
orem implies that the exponents i and j may be taken as residues modulo Icm/[3, 136] = 408 
(note that the calculations in Solution 2 imply that the multiplicative order of 2 modulo 17? 
is 136; this is not necessary for the construction of single examples as in that solution, but it 
is needed here in order to define the shape of canonical solutions). With these restrictions, 
we will show that there are exactly 51 unique solutions, which are listed at the end. 

In fact, all solutions lie in two simple families. Indeed, any (i, j) such that 2° + 2/ = —1 
(mod 2023) also projects to a solution modulo 17. The order of 2 modulo 17 is 8, and there 
are two solutions: (3,3) and (5,0). We now must check whether these lift to solutions modulo 
17°. 

In the first case, such a lift has the form 


peers  288+3 =] (mod 17°) 
for some integers 0 < a, 8 < 16. The Binomial Theorem implies that 


28% = (15-17 +1)%°=1+a-15-17 (mod 17°). 
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Plugging this in, the above equation reduces to 
(8(a+ B)-154+1)-17=0 (mod 17°) <=> 8(a4+)-15+1=0 (mod 17), 


which reduces to a+ 3+1=0 (mod 17). Since a, 3 play a symmetric role, the possible sets 
of solutions modulo 17? are spanned by choosing 0 < a < 8 and 6 = 16—a. To recover 
solutions modulo 2023, multiples of 136 must be added to (8a + 3,86 + 3) until the values 
are congruent to {1,2} (mod 3). For a < 8 there are two distinct ways of doing this, but 
for a = 8 there is only one, making for 17 total solutions. (In fact, when a = 8 we recover 
the congruence from Solution 2, namely 267 + 267 + 20 = 0 (mod 172), which lifts uniquely 
o (203, 67)). 

For the case (5,0), the calculations are similar. We find that 2°¢+5 +289 = —1 (mod 17?) 
has solutions parameterized by 8 = 1+2a mod 17. Now each such pair lifts to two solutions 
modulo 2023, since a and ĝ are no longer symmetric, so this gives a total of 34 solutions. 
For example, œ = 0, 8 = 1 corresponds to 2° + 28 + 2° = 1 (mod 177) from Solution 1, and 
a = 8,8 = 0 corresponds to 2° + 2° + 2° from Solution 2. 

The set of all reduced (i,j) is as follows: 


(77, 16), (85, 32), (101, 64), (109, 80), (125, 112), (133, 128), (139, 131), (155, 115), 
(160, 149), (163, 107), (176, 157), (179, 91), (187, 83), (203, 67), (208, 173), (211, 59), 
(224, 181), (227, 43), (235, 35), (251, 19), (256, 197), (259, 11), (272, 205), (277, 8), 
(280, 5), (293, 40), (296, 13), (301, 56), (304, 221), (317, 88), (320, 229), (325, 104), 
(328, 29), (341, 136), (344, 37), (347, 331), (349, 152), (352, 245), (355, 323), (365, 184), 
(368, 253), (371, 307), (373, 200), (376, 53), (379, 299), (389, 232), (392, 61), (395, 283), 
(397, 248), (400, 269), (403, 275). 
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B3. A sequence y1, Y2,..., Yk of real numbers is called zigzag if k = 1, or if yo — y1, Y3 — 


Y2,..-, Yk — Yk—ı are nonzero and alternate in sign. Let X1, X2,..., Xn be chosen indepen- 
dently from the uniform distribution on [0,1]. Let a(X1, X2,..., Xn) be the largest value of k 
for which there exists an increasing sequence of integers 71, 72,..., ip such that Xi, Xi,,... Xi, 


is zigzag. Find the expected value of a(X1, X2,..., Xn) for n > 2. 
Answer: (2n + 2)/3 for n > 2. 


Solution 1: We begin by noting that with probability one, the X; are all distinct. Indeed, the 
event in which there are at least two identical values defines a finite collection of hyperplanes 
in [0, 1)”, which has measure zero. Furthermore, since we consider only relative order , we can 
translate the problem to computing a(w), where w = w 1... Wn is a permutation of [1,...,7] 
chosen uniformly at random. 

Let u(w) be the length of maximal zigzag sequence starting with a descent, and let d(w) 
be the maximal zigzag sequence starting with an ascent. Then a(w) = max{u(w),d(w)} 
and moreover, for n > 2 we see that u(w) = d(w) +1, as we can obtain one from the 
other by adding/removing an initial element. By symmetry we see that u(w) and d(w) are 
identically distributed and exactly half of the time u(w) will be the larger. Thus, by linearity 
of expectation, we have 


Bfa(w)] = Elu(w)] + [Ell] = Blu(w)] + $ = Ela(w)] + 5. 

We claim that a maximal zigzag subsequence can always be chosen containing the value 
n. If a zigzag subsequence has w, = n where ij < r < 2;41, then we can reassign ij = r if 
Wi; > Wizz, OF 41 = 7 if wi; < Wi,,,, resulting either way in a zigzag subsequence of the 
same length containing n. If r < 71, we can reassign 7; = r if wi > Wi, or prepend wr to Wi 
if wi, < Wi, resulting an a zigzag subsequence of the same length or longer. The case r > ig, 
where k is the length of the subsequence, is similar. 

Let Sn denote the permutations of [1,..., n], and define 


fn := Eld(w)|w € Sn], gn := Ela(w)|w € Snl. 


By the above reasoning we have gn = fn + 5: 

Let w € Sni1. By the above reasoning, we can assume that the element n + 1 is part 
of a maximal zigzag sequence wj;,---wi,. Suppose that wj41 = n + 1 and ir = j + 1. Then 
Wi, < Wi, >+- is a maximal UD alternating subsequence (i.e. starting with an ascent) 
in the reverse permutation w;...w1, and Wip, < Wi,,, > ++: is a maximal UD alternating 
subsequence in wj+2...Wn+41. Thus we see that, choosing j with probability a we have 

1 
n+1 


Gn4+1 = 


(fj + fn—j + 1). 
j=0 


NIW 


We have the following boundary values, which are verified directly: fọ = 0, fi = 1, fo 
(d(12) = 2 and d(21) = 1). 
Then 


2fo+2fi +2% f| +1 


24 siete 1 
Gn+1 = Jn+1 5 peed 2 
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Rewriting this recursion for f we get 


n 


2 1 
— 1 - ae 
fn+1 nA F X Íj t3 


j=2 


Note that this implies a relation of fo +---+ fn via fn+1, which can be then reiterated 
for n — 1 and substituted 


HE (fost — 5) <1 fate ha t fn ES (h 5) = 5 (hm 5) + fin 


2 2 2 2 
Then 
(n+ Dfa = (M4 Diet E= 41) (on - 5) = (n+ 2) (o.- 5] + 
so 
(n+ Dgn = (n+ 2)gn => SEE = SR = =8=2 


which completes the proof. 


Solution 2: The event that two of the X,;’s are equal has probability zero, so in the argument 
below we assume that no two are equal. For 1 < j < n, call j a “turning point” if Xj — Xj~1 
and Xj; — X; have opposite signs, and let T be the total number of turning points. 

We claim that a(X1,...,Xn) =T +2. For 1 <j <n, let Dj = Xj41 — Xj. Suppose that 
Xj,,...,X4, is zigzag. For each 1 < m < k, we have Xj,,., — Xim = Dim + Dim+1 +: + 
Dinyi—1- Choose £m with im < lm < im+1 so that De„ has the same sign as Xj,,,, — Xi 
Then the sequence De, Dy,,..., De,_, alternates sign. Thus, there are at least k — 2 changes 
of sign in the sequence D1, D2,..., Dn-1. Each such change of sign is a turning point, so 
k-2<T,andk<T+2. To see that k = T + 2 is possible, let j1, jo,..., jr be the turning 
points. Then D1, Do,..., Dn-1 changes sign only between D,,, and D;,,-1 form = 1,2,...,T. 
It follows that the sequence 


m* 


Xj — X1, Xj, — Xis Xir — Xiri Xn = Xjr 


= Di +: + Dag Dji +++ + Diaa, Dira +t + Dir-1, Dir +++ + Dn- 


alternates sign, so X1,.Xj,,...,Xj,-,Xn is a zigzag subsequence with length T + 2. This 
verifies our claim. 

For 1 < j < n, the probability that j is a turning point is 2/3, since 4 of the 6 equally likely 
orderings of X;j—1, Xj, Xj+1 yield a turning point. Therefore, E[a(X)] = (n — 2)2/3 + 2 = 
(2n + 2)/3. 


Solution 3: Let Zj = a(Xj,...,X;) — a(X1,...,X;-1) for j > 3. Then for n > 2, we can 
write a(Xy,...,Xn) = a(X 1, X2) + Z3 +--+ + Zn, so it suffices to determine 


Bla(X1, X2) +Z +: + Zn] = Eļa(Xı, X2)| + E|Z3] ig E[|Zn]. 


Notice that Eļa(X1, X2)] = 2 since with probability one either Xı > Xə or Xı < Xə, and in 
either case X1, X2 is zigzag. We claim that E[Z;] = 2/3 for j > 3, from which the answer 
2 + 2(n — 2)/3 = (2n + 2)/3 follows immediately. 
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To verify the claim, we first excluded the possibility that two of the X;’s are equal, 
which has probability zero. Next, notice that if Xj;,,...,Xi, is zigzag and 7 > i, then 
either X;,,...,Xi,_,,Xi or Xi,,...,Xi,, Xi is zigzag, according (if k > 2) to whether or 
not X;, is between X;, , and X;. Then if k = a(Xj,...,Xj-1) and Xj,,..., Xi, is zigzag, 
Xj,,...,Xi,_,,Xj-1 must also be zigzag (since otherwise there would be a longer zigzag 
subsequence). Further, X;_; — X;-2 must have the same sign as X;_; — X;,_,; otherwise, 
Xj-2— Xi, = (Xj-2— Xj-1)+ (Xj-1— Xi,_,) would have the same sign as Xj—2—X;—ı and 


k-1 k-1 
Xj—1 — Xi,_,, hence the opposite sign as X;_; — X;_2, making a longer zigzag subsequence 
A eas jie Xj—2, Xj-1. Thus, Zj = 0 if Xj = Xj-1 has the same sign as Xj-1 E Xj-2, 


and Z; = 1 if they have opposite signs. In other words, Z; = 0 if X;_1 is between Xj_2 and 
Xj, and Z; = 1 otherwise. Since X;_1 is between Xj_2 and X; for 2 of the 6 equally likely 
orderings of X;j—2, X;~-1, X;, the probability that Z; = 1 is 4/6 = 2/3, and E[Z,] = 2/3 as 
claimed. 
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B4. For a nonnegative integer n and a strictly increasing sequence of real numbers tg, t1,..., tn, 
let f(t) be the corresponding real-valued function defined for t > to by the following proper- 
ties: 


(a) f(t) is continuous for t > to, and is twice differentiable for all t > to other than t),..., tn; 


(b) f(to) = 1/2; 


(c) lim, ,¢+ FA =0fr0<k<n; 


(d) For 0 < k < n-— 1, we have f”(t)=k+1 when t, < t < tk+1, and f”(t) =n +1 when 
t > th. 


Considering all choices of n and to,t1,...,t, such that tk > tk-1 +1 for 1 < k < n, what is 
the least possible value of T for which f (to + T) = 2023? 


Answer: 29 


Solution 1: Let T be the value for which f(tp + T) = 2023, and assume without loss of 
generality that tn < to + T, since greater values do not affect f (to +T). Let tn+1 = to + T. 
Notice that for each 1 < k < n+1 we have f(t) = f(t,-1) +k(t—tp_1)?/2 when ty_1 < t < tp. 
Let Tk = tk — tk—1; then Tk > 1 for 1 < k < n and T41 > 0. Let m = n + 1. The goal is to 
minimize T = 7, + T2 + -+-+ Tm subject to the constraint 


C(t) = So kT = 2(f (to + T) — f(to)) = 4045. 
k=1 


The space to be minimized over consists of all m > 1 and all real m-tuples (T1,...,Tm) 
with 7 > 1 — dy, where mk is the Kronecker delta. This space can be made topologically 
connected with the identification (71,...,Tm-—1,0) = (T1, - --, Tm—-1) for m > 2. The subset of 
this space that satisfies the constraint is bounded and hence compact, because the constraint 
excludes values of m for which m(m — 1)/2 > 4045 and values of T greater than /4045. 

Let (71,72,°:* ,7m) lie in the constrained space, so that C(r) = 4045 and Tk > 1 — bmp. 
If there is more than one value of k for which Tk > 1 — dmx, let two of these values be kı 
and k2, and assume without loss of generality that ki7,, > k2Tķọ. Choose £e > 0 sufficiently 
small that Tk, — € > 1—dms., and consider the m-tuple for which (Tk, Tka) is replaced by 
(Tki +E, Tka — €), and the other 7 are unchanged. Then T is unchanged, while the sum 
> 7, kr? in the constraint changes by 


2(kiTk, — k2Tko )E + (ky + ka)e?, 


which is (strictly) positive by our earlier assumption. Thus, the sum is now greater than 
4045. Next, reduce Tk, + € until the constraint is satisfied again, and notice that the value 
that meets this condition is strictly between Tk, and Tk, +€ (since C(r) is increasing in each 
Tk), so it satisfies the same lower bound as Tk. The resulting m-tuple has a strictly lower 
value of T, so the original m-tuple could not have minimized T. Thus, the (constrained) 
minimum of T cannot be achieved with more than one value of k for which 7 > 1 — dyn. 
Next, suppose that T is minimized with Tọ > 1 for some k < m; then by the argument 
above, Tm = 0. In that case, reduce m by 1, and since Tm—1 > 0, the above argument implies 
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that T is not at a minimum for the new value of m, and hence not a global minimum. By the 
compactness described above, there must be a global minimum value for T, and by what we 
have argued so far, this minimum must satisfy 7, = T2 =--- = Tm—-1 = 1. Thus, the problem 
is reduced to minimizing 


m—-1 


T(m) =m—1+4+ (1045-5 6) m= m1 + VERT HT 


k=1 


over all 1 < m < M, where M is the greatest integer for which the square root is well-defined. 
The fact that M > 20 suffices for the arguments below. 
If m < M, then T(m +1) > T(m) is equivalent to 


1+ /4045/(m + 1) — m/2 > ./4045/m — (m — 1)/2. 


Squaring both sides (which we observe are nonnegative) and simplifying, this is equivalent to 


2,/4045/(m + 1) — m/2 > 4045/[m(m + 1)] — 1/2. 


Notice that the right side is nonnegative when the left side is well-defined (when m < M). 
We square and simplify again to conclude that T(m + 1) > T(m) is equivalent to 


4m > 4045/[m(m + 1)] — 1/2 => m(m + 1)(m + 1/8) > 1011.25. 


Thus, T(m + 1) > T(m) for 10 < m < M. 

Similarly, T(m + 1) < T(m) is equivalent to m(m + 1)(m + 1/8) < 1011.25, which is true 
for 1 < m < 9. We conclude that T(10) = 9 + v 404.5 — 4.5 = 29 is the minimum possible 
value of T. 


Solution 2: Following the same notation and initial steps as the previous solution, the goal 
is to minimize T = 71 + T2 +--+- +7 41 Subject to the constraint 


n+l 
C(r) = So kr = 4045, 
k=1 
where n > 0 and Tk > 1 for k = 1,...,n and ™m41 > 0. 
We now make the following substitution: Let £k = ke — k > 0 for k = 1,...,n and 


Chi = (n+ 1ra- We have x; > 0 and 
1 
In41 = 4045 — e J-a- mao, 


Then the goal is to minimize 


n n+1 
lop. +k 4045 £i — - x 
Pi E E - y aie r 
k=1 


over the simplex x; > 0 and 71 +---+ £n < 4045 — (20) and over n > 0 for which this 
simplex is non-empty. (We include n = 0, for which the simplex is undefined but the original 
constraint requires Tı = 4045; in this case, T = 4045, so we let Fo = v 4045.) 
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We now analyze the function gk(£k) = Fn(£1,..., £n) for fixed values of zj, j Ak. We 
have 


, 1 1 
ATK) = 
2y kler +k)  2,/(n + 1) (4045 — ($1) — z1 — +- — £n) 
i 1 1 
9 (Lk) = <0, 
À AVk (a +k)3/? dyar a= aa 


so gk is a concave function and the minimum is achieved at the boundary. Thus we either 
have zk = 0 or 41 +: + £n = 4045 — (Oe) for every k. So the minimum is achieved for 
some (21,...,2n) either on the hyperplane Hn := {£1 +---+ £n = 4045 — Coot or else we 


must have x, = 0 for all k so (£1,..., £n) = (0,...,0). 
In the first case we have x, = 4045 GS) £1 — +++ — Lpn_1, SO 
n—1 n+1 
fx, +k 4045 — Ti == 

TTE S T Fi ( 2 ) = 4 = Ppi lher): 

k=1 
If (£1,...,&n—1) minimizes Fp—1 and (z1,...,£n—-1) € Hn-1, then we can again reduce to the 
case of one fewer variable. Recursively, we can eventually reach a minimum in the second 
case where (£1,...,£Ẹ) = (0,...,0), which is vacuously true for N = 0 in case we get that 
far. 


Having reduced to the second case, we now need only minimize the function 


4045 — (Nt?) 4045 
t(N) := 0,...,0)=N4 
(N) = fin ) we EUNET 


over the integers 0 < N < 89 (since the square root is undefined for N > 90). To motivate 
the next steps, crude estimates show that t(1) > 40, t(3) < 40, t(10) < 40, a t(30) > 40. 


This suggests that the minimum occurs when Wet is considerably larger than x To get 
a rough idea where the minimum occurs, replace N for the moment with a real variable y, 


a . From its derivative, the latter function is 


and consider the approximation t(y) ~ y + 
minimized when y = —1 + (4045/4)!/3 ~ 9. 
Consider then t(y) — t(9) = t(y) — 29; we will find values of y > 0 for which this is 


nonpositive. Thus, we are solving for 


4045 4045 y 
n Z <29- y = y < 29 and —— — Z < (29 — 
T a TA ee 


The last inequality is equivalent to (using the fact that we have equality for y = 9 to factor) 


(y — 29)? (y + 1) — 4045 4 wy D L 1 ory — 9)(2y? — 95y + 712) 
R li 9) ( 95 + =) ( 95 — =) si 


This inequality is true only when y > Sty 3329 3329 ~ 29 or when y is (not strictly) between 9 
and 95=73329 y3329. The latter number is ay between 9 and 10 because 55? = 3025 < 3329 < 


3481 = 59°. Thus, t(y) — t(9) < 0 only when 9 < y < = < 10. Then ¢(N) > t(9) for 
integers N, and t(9) = 29 is the desired minimum. 
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B5. Determine which positive integers n have the following property: For all integers m that 
are relatively prime to n, there exists a permutation 7: {1,2,...,n} —> {1,2,...,n} such 
that (a(k)) = mk (mod n) for all k € {1,2,...,n}. 


Answer: Those n of the form 47 + 2 where 7 is a nonnegative integer, and n = 1. 


Solution: Suppose that T E€ Sn, where Sn is the group of permutations on n elements. We 
say that a permutation 7 € Sn is a square root of T if n? = r (we will shortly see that in 
general there is not a unique square root, if one exists). If m is an integer coprime to n, let 
Tm © Sn denote the permutation defined by multiplication by m modulo n, in other words, 
Tm(k) =m-k mod n for 1 < k< n. The problem statement asks for a classification of the 
n such that all Tm have a square root in Sy. 

Now consider the standard decomposition of 7 into disjoint cycles. For a positive integer 
i such that i < n, let f;(7) denote the number of distinct cycles of length i. 


Lemma. Suppose that T E€ Sn. Then T has a square root if and only if foi(7) is even for all 
a< 8, 
=9 


Proof. If £ is a positive integer, consider the permutation in S, defined by a cycle of length 
L, say y = (x£1£2--- £e). When squared, the resulting permutation either preserves the cycle 
length, or splits it into two cycles of half the length, depending on the parity of £: 
zi _ (x1 23 +++ Lok410204 °° £93) if @=2A+1, (3) 

(£1£3-:: £2)—1)(£2£4 <- £2) if C= 2A. 


Applying this fact to the cycle decomposition of a permutation m € Sn implies that (noting 
that all of the resulting cycles remain disjoint) 


forgi(m) = forgi(m) + 2faa+2(7), 
faln’) = 2 faa (T). 


Here we set f;(7) = 0 if i >n. The forward direction of the claim follows, as f(r?) is even 
for all even i. 
For the reverse direction, suppose that 7 is a permutation such that fa;(7) is even for all 
i. We now construct a square root 7 as a list of cycles, essentially using (3) in reverse. For 
each cycle in 7T of odd length, say y = (£1£2--- £2)+1), append the following length 2A + 1 
cycle to 7: 
(1%) 49%2%y43%3 ++ T2)+1TA+1)- 


By assumption, there are an even number of cycles of even length in 7, so they may be 
grouped (arbitrarily) in pairs. For each pair of the form (x12: £2)), (Y1 **-Y2)), append 
the following cycle of length 4A to a: 


(£1Y1£2Y2 is T2AY2A)- 


Now it is immediate that 7? = 7, since by construction they share the same cycle decompo- 
sition. 
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This result also shows that even parity of 7 is a necessary condition for the existence of 
a square root (though not sufficient, as for example, T = (12)(3456) is an even permutation 
without a square root). 


Corollary. If 7 is an odd permutation, then it has no square root. 


Returning to the problem at hand, note that the cycles of Tm € S, are of the form 
(m) -x := (x mg -mT c), 


where x is a positive integer, and @ is the minimum positive exponent such that m‘x = x 
(mod n). This is equivalent to the minimal £ such that mf = 1 (mod AGA 

If n = 4k + 2, then consider an arbitrary value of m, which must be odd. In this case, 
we will show that fi(Tm) is even for all i, and therefore Tm has a square root by the Lemma. 
Specifically, we claim that fi(Tm) = 20;(Tm), where 0;(Tm) denotes the number of cycles of 
length 2 such that all elements are odd. 

To verify this claim, suppose that x is odd. Then (m) -x consists of only odd integers 
(and any such cycle is generated by some odd x), and maps bijectively to (m) - 2x (a cycle 
with only even integers), since x ++ 2x (mod n) bijectively maps the odd residues modulo n 
to the even residues (noting that n is divisible by 2 but no higher power of 2). 

If n = 1, then the property trivially holds because all integers are congruent to each other 
modulo 1, and because there does exist a (trivial) permutation 7: {1} > {1}. 

For all other n, we will provide a value m = mp such that Tm„ does not have a square 
root in Sn. If n = 4k, let m, = —1. Then fi(t) = 2 (the fixed points being z = 0 and 
2k) and f2(Tm) = 2k — 1, so the conclusion follows by the Lemma. A similar argument also 
shows that 7_; does not have a square root if n = 4k + 3. However, this does not work for 
n = 4k +1, and the case that n is odd can be treated in a unified manner as follows. 

Suppose that n > 1 is odd, with prime factorization n = p{'---p", where the p; are 
distinct odd primes, and e; are positive integers. Let mı be a primitive root of the multi- 
plicative group modulo pî! (cf. Euler, Lagrange, Legendre, Gauss for the existence of such 
a root). Then mı is also a primitive root modulo p{ for 1 < e < eı — 1. Now (using the 
Chinese Remainder Theorem) set m to be the residue modulo n that satisfies 


m=m, (mod pt), m=1 (mod p;’) for 2 < j <r. 


Then all non-fixed-point cycles of Tm have lengths of the form py (pı —1) for some 1 <e<e. 
Specifically, if x is a multiple of p}', then (a) -x has length 1, and Pra if pj is the largest 
power of pı dividing x (where 0 < s < e; — 1), then the length is pft ° (pı — 1). 

Thus by considering the cycles formed by all x’s that are not multiples of pı (this is the 


case s = 0, but in fact any s < e; works similarly), we find that 


Since this is odd, the Lemma implies that Tm is not a square. 
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Remark. The Corollary can also be proven directly (almost immediately using the fact that 
the sign map is a homomorphism). This gives an alternative criterion for showing that T—1 
does not have a square root when n = 0,3 (mod 4). 

There are many possible constructions for permutation square roots in the case n = 4k+4+2. 
The proof above implicitly defines a square root of Tm by “zipping” together all cycles of the 
form (£, mg, m?zx,...) and (22,m-2x,m?-2zx,...). One alternative is to instead pair cycles 
of the form (x,ma,m?z,...) and ((a+n’),m(a+n’),m?(x+n’),...), where n’ = 2k+1 = 3. 
To be more precise, it is straightforward to show that the following permutation is also a 
square root of Tm: 


n(x) := 


x+n' (mod n) if x is odd, 
mr+n (mod n) if x is even. 


28 


B6. Let n be a positive integer. For i and j in {1,2,...,n}, let s(i, j) be the number of 
pairs (a,b) of nonnegative integers satisfying ai + bj = n. Let S be the n-by-n matrix whose 


(i, j)-entry is s(i, 7). 


6 3 2 2 2 
3 0 1 0 1 
For example, when n = 5, we have S= | 2 1 0 0 1]. 
2000 J 
i 1 1 1 2 


Compute the determinant of S. 


Answer: (—1)®+t12m where m is the least integer greater than or equal to n/2 


Solution 1: Let d(i,j) = 1 if i | j and 0 otherwise, and let D = Gea) 4 o) be the 
corresponding n x (n +1) matrix whose rows are indexed by 1,...,m and whose columns are 
indexed by 0,1,...,n. For example, for n = 8, we have 


j=0 j=1 j=2 j=3 j=4 j=5 j=6 j=7 j=8 


i=l 1 1 J: 1 1 1 1 1 ik 
1=2 1 0 1 0 1 0 1 0 1 
i=3 1 0 0 1 0 0 1 0 0 
D=|i=4 1 0 0 0 1 0 0 0 1 
i=5 1 0 0 0 0 1 0 0 0 
i=6 1 0 0 0 0 0 1 0 0 
i=7 1 0 0 0 0 0 0 1 0 
i=8 1 0 0 0 0 0 0 0 1 
All possible solutions (a,b) to ai+bj = n in nonnegative integers can be indexed by those 


0 < k <n such that i | k and j | n — k, so 


n 


s(i, j) = X di, k)d(j,n — k). 
k=0 
Thus S = DMs. where D” is the matrix D with the rows reversed, i.e., D}; = Din—j;- 


UJ 
The Cauchy—Binet formula implies that if M is ann x(n + 1)-matrix and N is an (n + 1) x 


n-matrix, then 
det (MN) = Y” det (Mz) - det (n) l 
k=0 


where M; denotes the matrix M with the k-th column removed, and N k denotes the matrix 
N with the k-th row removed. 
In the present case we have M = D and N = (D")", and can further reduce 


det Cp = det ((D")z) = (—1)"/! det D~. 


The first equality is due to the fact that determinants are preserved by transposition, and 
the sign arises from writing the row reversals as |n/2| column swaps. In all, we therefore 
have 
n 
2 
det S = X` det Dz(—1)'"/?! det D~. (4) 
k=0 
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This sum can be further simplified by evaluating det D> for large k. 
Lemma. (a) Ifk > n/2, then det Dz = (—1)*". 
(b) If n is even, then det D~; = 0. 


/ 
Proof. For k > n/2 we can compute det Dz by cofactor expansion along row k. For such a k 
we have d(k, j) = 0 unless j = 0 or k, and on removing column k, row k of De is [1,0,..., 0]. 


Thus det Dz = (—1)*—!, since after removing row k and columns 0 and k of D, the remaining 
matrix is upper triangular with 1s on the diagonal. 

If n is even, say n = 2m, then row m of Dy is [1,0,...,0,1] (since d(m, j) = 1 only when 
j =0,m,2m, and j = m has been removed). Thus d(m, j) = d(2m, j) for j 4 m, and since 
row m and row n = 2m are then identical, det Dz = 0. 


Plugging in to (4), we can hence in all cases set m := [n/2] and write 


m—1 


det S = 2(-1)!"/7! X` det Dz det D = 2(-1)l/?! xt )"-*-1 det Dz 
k=0 
m-1 
= 2(-1)""! S°(-1)* det Dz, 
k=0 


where we have also used that n — m = |n/2|. The proof is complete on evaluating this sum. 


Claim. We have 


T 


m— 


E det (D 


Z 
3 


2l 


Proof. Let D’ be the (n + 1) x (n + 1)-matrix obtained from D by attaching an extra row 
(1,1,...,1) at the top of the matrix (i.e., at row index i = 0). Then, the topmost two rows of 
the matrix D’ are equal, so that det (D’) = 0. On the other hand, expanding the determinant 
of D' along the topmost row, and using part (b) of the Lemma to eliminate the middle term 
of the sum when n is even, we obtain 


n m-—l1 n 
det (D') = 5 (—1)* det (Dz) = X- (-1)" det (Dg) + XO (-1)* det (Dz) 
k=0 k=0 k=n—m-+1 
m—-1 n m—l1 
= X (-1)* det (D3) SO (-1)*- (-1)*"1 = X` (-1)* det (Dg) - 
k=0 k=n—m-+1 k=0 


The second line follows by part (a) of the Lemma. Recalling that det(D’) = 0, the proof is 
complete. 


Solution 2: As in Solution 1, the (i, j)—th entry of S is the convolution of indicator series 
for multiples of i and j. In particular, this gives the factorization S = B-Rev(B)", where 
Rev(B) denotes the matrix formed by reversing each row of B, and B consists of the rows 
Ri,..., Rn, where 

R; = (10---010---01---) =(1( | j))j=0 


i—l1 i-1 
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The key idea in this proof is to apply row operations directly on B. In particular, if M 
is an n x n matrix with det(M) = 1, then 


det (MB -Rev(B)'M*) = det(S). 


Furthermore, Rev(B)’M7 = Rev(MB)’, since under left-multiplication M acts as row 
operations on B, which are commutative with row reversal. 

The preceding claims can also be justified more explicitly by noting that RevB = M - J, 
where 


0 0 1 

0 1 0 
J := ; 

1 0 0 


For example, commutativity follows since M - Rev(B) = M - B - J = Rev(M B). 
The rows R; form a full rank system, since B excluding the first column is upper triangular. 
We can therefore reduce B to the canonical form 


ay 100. 0 0 
a 0 1 0 0 0 
MB = : : 
anı 0 0 0 `. 10 
an 000 >. 0 1 


by a sequence of row operations that acts on the rows from bottom to top, with each operation 
subtracting a row from below the row being acted on. It follows that a; can be defined 
recursively for i = n,n — 1,...,1 by the formula 


ai = 1 — ag — @3i— + -- 


where for convenience we set the initial conditions a; = 0 for i > n. Equivalently, we can set 
additional initial conditions aņn-m+1;---,@n = 1, where m = [4], and define a; recursively 
for i =n—m,...,1. 

Assume now that n is odd, so that n = 2m — 1; it will turn out that n = 2m reduces to 
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this case almost immediately. Then n —m + 1 = m, and 


MSM? = MB - Rev(B)' MT 


ai 1 0 0- 0 0 i : 1 , 
a2 0 1 0 -:--- 0 O 
2 fio . a 0 0 0 1 0 0 
1 0 > 1 0 
0 1 0 
1 0 0 0 0 1 aks o 
ay, ag ümar. l 1 
1 ay 
1 Am—-1 
= 1 1 
1 1 
Gp. Sek oaa A wee FE 2 


The determinant of this matrix can now be evaluated by reducing the final row to zeroes in 
all except the final column, and then expanding along that row. This gives 


det(S) = det( MSMT) = (-1)™1! (2 — 2(a1 +++» + am-1)). 


Lemma 4. With a; defined as above, 
ay +t + am-1 = 


Proof. Plugging in the definition of a, and initial 


—(m-—1). 


values for am, ..., an, 


ay +: + am- = (1 — a2 — a3 


an) + a2 +: + am-1 


=I aman =1-—™M. 


This completes the proof for n = 2m — 1. 


Finally, in the case that n = 2m, the recursion gives am = 0. Thus, the mth row of MB 
has only one nonzero entry, a 1 in the first column, and is orthogonal to all other rows of 
MB. It follows that the mth row in MSMT has a 1 in the mth column, and is 0 elsewhere, 
so this row has no effect on its determinant. Further, the conclusion in Lemma 4 also still 
holds, since now am = 0 implies that am + a@m41 +++: + an =m. 


Solution 3: As in Solution 2, let m be the least integer greater than or equal to n/2, so that 


either n = 2m or n = 2m — 1. We will define (n 4 
that 
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—2m O 


+ 1) x (n+1) matrices C and D and show 


s | 


where O denotes a row of n zeros. It will then follow that det S = (det C)(det D)/(—2m). 
We will write the entries of C and D as cij and dij where0 <i<nand0<j <n. 

For 1 <i <n and0 < k <n, let cip = 1 if k is a multiple of i, and cj; = 0 otherwise. For 
O<k<nand1<j <n, let dg; =1ifn—k is a multiple of j, and dy; = 0 otherwise. Then 
each allowed solution of ai + bj = n corresponds to a case where Cik = dkj = 1, with k = ai. 
Thus, s(i, j) = Th Cikdkj, verifying the S block in the equation above for CD. 

Next, we will choose column 0 of D (corresponding to j = 0) to be orthogonal to rows 1 
to n of C, which will satisfy the OT block of the equation for CD. Such a column must exist 
because these n rows can’t span (n + 1)-dimensional space. Notice that the n x n matrix 
{Cik}1<i,k<n is upper triangular, with all ones on its diagonal, so the rows of this reduced 
matrix are linearly independent. Thus, for the desired orthogonality, dog must be nonzero, 
and since rows 1 to n of C span an n-dimensional space, choosing doo = 1 uniquely determines 
column 0 of D. For n—m +1 < i < n, we have 27 > n, so cig = 1 if and only if k = 0 or 
k = i. Thus, the desired orthogonality requires that diọ = —1 for these values of i. Also, if 
n is even, then row m of C is the same as row n of C except that Cmm = 1 while Cnm = 0. 
Thus, the desired orthogonality requires that dmo = 0 when n is even. Then, whether n is 
even or odd, we have X` }—m dko = —m. For column 0 of D to be orthogonal to row 1 of C, 
all of whose entries are 1, we then must have ies dko = m. We will not need to determine 
the individual values of do, for 1 <<k<m-—1. 

Now let cox = dn—zo for 0 < k < n. Notice that this makes row 0 of C (which is the 
reverse of column 0 of D) orthogonal to columns 1 to n of D (which are the reverses of rows 
1 to n of C), satisfying the O block of the equation for CD. To complete the verification of 
this equation, the upper left entry of CD is 


n n m-1 n 
S cokdro = X dn—odko = X (—1)dro+ XO da—zo(-1) = -m—m = —2m. 
k=0 k=0 k=0 k=n—m+1 


Next, to determine det C', replace column 0 of C with column 0 of CD. Since the latter 
column is the sum over 0 < k < n of dko times column k of C, and dog = 1, this replacement 
is a column operation on C that this does not change its determinant. The resulting matrix is 
upper triangular, and its diagonal consists of —2m followed by n ones. Thus, det C = —2m. 

Finally, since D can be obtained by reversing the columns of CT, which amounts to swap- 
ping m pairs of columns, det D = (—1)™det CT = (—1)™ det C = (—1)™+!2m. Therefore, 
det S = —2m(—1)™*12m/(—2m) = (-1)'™*12m. 
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